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ABSTRACT 


A  method  for  determining  the  complex  bulk  modulus  B  and  the 
complex  mass  density  p  of  a  fluid  contained  in  a  rigid  porous  solid  is 
investigated.  The  fluid  filled  solid  is  contained  within  a  small 
cylindrical  cavity  capped  on  each  end  by  identical  transducers.  This 
method,  suitable  for  an  acoustically  "soft"  transducer,  is  based  upon 
electroacoustic  network  theory.  B  and  p  are  obtained  from  the  input 
electrical  impedance  of  each  transducer  when  they  are  wired  in 
parallel  "in-phase"  and  "out-of-phase",  respectively.  A  description  of 
the  apparatus  and  preliminary  experimental  results  are  presented. 


f  OTIC 

corr 

cry 


j  Accession  For 

V 

j  NT  18  GRA4I 

fflT 

;  OTIC  T:'.F 

n 

!  tli  X.  c  ed 

n 

•J  *  •  '  t  C t  j  OH _ 

i  Dl  .-jtrlbufcloe/ 

I  AvJSiiftbiXlt;'  ’-liC'i 

•  Du>C  '  ?.ij\  •: ii.i 


ill 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

n.  THEORY . 3 

A.  INTRODUCTION . 3 

B.  ACOUSTIC  IMPEDANCE  IN  SAMPLE  CAVITY . 3 

1.  Introduction . 3 

2.  Blocked  Acoustic  Impedance . 5 

3.  Transfer  Impedance . 7 

4.  Long  Wavelength  Approximation . 8 

5.  Identical  Transducers . 10 

6.  Fluid  Filled  Porous  Solid . 1 1 

C.  TRANSDUCERS . 11 

D.  COUPLER-DRIVER  SYSTEM . 14 

III.  TRANSDUCER  CHARACTERIZATION  AND  COUPLER . 1 7 

DESIGN 

A.  BACKGROUND . 17 

B.  TRANSDUCER  SELECTION  AND  PARAMETER . 17 

CHARACTERIZATION 

1 .  Matched  Transducer  Identification . 17 

2.  Source  Strength . 1 9 

3.  Transducer  Parameter  Determination . 21 

a.  Electrical  Parameters . 21 

b.  Driver  Stiffness . 23 

c.  Moving  Mass . 28 


iv 


d.  Transduction  Coefficient . 30 

e.  Mechanical  Resistance . 33 

f.  Cavity  Stiffness . 33 

g.  Modeled  Performance . 36 

C.  COUPLER  DESIGN  AND  CONSTRUCTION . 39 

1.  Coupler  Construction . 39 

2.  Fluid  Management  System . 41 

IV.  PRELIMINARY  EXPERIMENTS  AND  RESULTS . 44 

A.  INTRODUCTION . 44 

B.  WATER-FILLED  TUBE . 44 

1.  Description  of  Tube . 44 

2.  Resonance  Experiment . 44 

3.  Cross-Sectional  Area  Correction  Factor . 45 

4.  Broad  Band  Experiment . 48 

C.  COUPLER  EXPERIMENTS . 50 

1.  Repeat  of  Water  Filled  Tube  Experiment  With  Half....  50 

of  Coupler 

a.  Resonance  Shift . 50 

b.  Area  Correction  Factor . 51 

c.  Broad  Band  Experiments . 51 

2.  Experiments  With  Complete  Coupler . 5 1 

a.  Air  filled . 51 

b.  Water  Filled . 54 

V.  CONCLUSIONS  AND  RECOMMENDATIONS . 57 


v 


APPENDIX . 58 

REFERENCES . 62 

INITIAL  DISTRIBUTION  LIST . 64 


vi 


ACKNOWLEDGEMENTS 


Several  persons  and  institutions  have  aided  me  by  providing 
support  and  assistance.  In  particular,  I  appreciate  the  support  of  the 
Naval  Research  Lab  Underwater  Sound  Reference  Detachment, 
Orlando.  I  thank  Professors  Wilson  and  Baker  for  their  patience  and 
guidance.  Finally,  1  treasured  the  encouragement  of  my  sons,  the 
understanding  of  my  wife,  and  I  thank  God  all  this  could  happen. 


I.  INTRODUCTION 

Geoacoustic  models  of  fluid -saturated  sediment  contain  several 
parameters  which  depend  on  the  microscopic  pore  geometry  of  the 
sediment  [Ref.  1,  Ref.  2,  Ref.  3].  These  parameters  are  very  difficult  to 
determine  from  the  usual  measurements  of  dilatational  sound  speed 
and  attenuation.  In  particular,  the  parameters  which  determine  what 
is  essentially  the  complex  effective  mass  density  (p)  of  the  fluid 
contained  in  the  pores  are  among  the  least  well  known  [Ref.  4,  Ref.  5]. 
For  this  reason  it  would  be  very  desirable  to  directly  measure  the 
complex  effective  fluid  mass  density  and  thereby  extract  these 
parameters  for  various  typical  microscopic  pore  geometries.  To  date, 
direct  measurements  of  p  have  been  conducted  with  helium  filled 
sintered  bronze  samples  [Ref.  6[. 

A  method  is  proposed  for  measuring  the  complex  effective  mass 
density  of  a  fluid  contained  in  a  rigid  porous  solid.  The  fluid-filled 
solid  is  contained  in  a  cylindrical  cavity  capped  on  both  ends  by 
identical  transducers.  The  method  investigated  is  an  impedance 
method  based  upon  electroacoustic  network  theory.  This  method  may 
also  be  used  to  measure  the  complex  bulk  modulus  (B),  and  thereby 
obtain  the  complex  speed  of  sound,  using  only  a  small  sample.  The 
impedance  method  is  suitable  for  use  with  acoustically  "soft" 
transducers  at  the  ends  of  the  sample  chamber.  Moving  coil,  high- 
fidelity  tweeters  were  used  in  this  investigation.  B  and  p  are  obtained 
from  the  input  electrical  impedance  of  the  transducers,  which  is 
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related  to  the  acoustic  impedance  presented  to  each  face,  when  they 
are  wired  in  parallel  "in-phase"  and  "out- of- phase",  respectively.  If  the 
two  transducers  are  identical  in  electrical,  mechanical,  and  acoustical 
characteristics,  data  reduction  is  remarkably  simplified. 

The  objectives  of  this  research  are  to  (1)  identify  and  characterize 
two  identical,  acoustically  "soft"  transducers;  (2)  develop  a  compact 
apparatus  to  measure  both  B  and  p;  and  (3)  validate  the  technique 
using  samples  for  which  the  desired  values  are  known  or  easily 
calculated. 
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II.  THEORY 

A.  INTRODUCTION 

In  this  chapter  the  network  equations  from  which  B  and  p  can  be 
extracted  are  developed  from  the  input  electrical  impedance  of  the 
transducers.  This  development  is  presented  in  three  major  sections. 
In  the  first  section  the  network  theory  associated  with  the  acoustics  of 
the  chamber  is  presented.  This  development  is  general  in  form  but  is 
simplified  by  use  of  the  long  wavelength  approximation.  The  specific 
case  of  matched  transducers  is  then  analyzed,  resulting  in  a  simple 
expression  for  the  acoustic  impedance  presented  to  the  face  of  each 
transducer. 

In  the  second  section  the  network  theory  for  the  transducers 
selected  for  the  experiment  is  presented.  The  relationships  between 
the  electrical  and  mechanical  parameters  are  developed.  These 
parameters  are  then  modeled  in  an  equivalent  electrical  circuit.  The 
third  and  final  section  consolidates  the  acoustic  network  of  the 
chamber  and  the  electroacoustic  network  of  the  drivers.  After  this 
consolidation  is  accomplished,  circuit  theory  is  then  used  to  derive  B 
and  p  from  the  measured  input  electrical  impedance. 

B.  ACOUSTIC  IMPEDANCE  IN  SAMPLE  CAVITY 

1.  Initial  Development 

The  goal  of  this  section  is  to  demonstrate  how  the  two  complex 
parameters  B  and  p  can  be  extracted  from  the  acoustic  impedance 
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presented  to  the  face  of  the  transducers  on  each  end  of  a  fluid-filled 
porous  sample  in  a  cavity.  A  simple  diagram  to  aid  in  explaining  this  is 
shown  in  Figure  1. 
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Figure  1 .  Simple  Chamber. 


In  this  figure,  there  are  two  transducers,  or  drivers,  facing  one 
another  with  the  sample  in  between.  The  cross-sectional  area  of  the 
drivers  are  assumed  to  be  equal  to  the  cross-sectional  area  of  the 
chamber.  The  two  port  network  representing  the  acoustic  properties 
of  this  chamber  is  shown  in  Figure  2. 


Figure  2.  Acoustic  Two  Port  Network. 


In  Figure  2,  all  the  parameters  listed  are,  in  general,  complex.  A 
complex  parameter  is  represented  with  a  tilda  (~)  over  the  symbol. 


The  symbols  are  defined  as: 

p  =  the  acoustic  pressure  at  the  face  of  the  transducer. 

U  =  the  volume  velocity  produced  by  the  transducer. 

ZQ  =  the  acoustic  impedance  seen  by  a  transducer  when 
the  other  transducer  is  held  fixed,  or  blocked. 

T=  the  transfer  impedance  between  the  two  transducers. 

The  appearance  of  Z0  and  T  in  each  branch  of  the  network  is  a 

consequence  of  the  reciprocal  property  of  the  acoustic  field  and  the 
symmetry  of  the  chamber.  Expressions  for  ZQ  and  T  for  can  be 

extracted  when  the  appropriate  boundary  conditions  and  wave 
functions  are  applied,  as  described  below.  The  subscripts  1  and  2  In 
Figure  2  associate  the  parameter  with  its  respective  driver. 

2.  Blocked  Acoustic  Impedance 

The  acoustic  pressure  and  volume  velocity  can  be  obtained  using 
the  concept  of  a  velocity  potential  (O).  If  the  velocity  potential  is 
known,  the  acoustic  pressure  (p)  and  particle  velocity  (u)  are  derived 
using  the  following  relationships  [Ref.  7:pp.  104-106] 

u=VO,  (2.1) 


dO 

p  =  -pG  —  .  where  p0  is  bulk  density. 

at 


(2.2) 


As  an  example,  for  a  time  harmonic,  sinusoidal  wave  in  one  dimension 
(  u  can  now  be  expressed  as  u  )  of  the  form  e^1  ^  ,  these 


expressions  reduce  to 

u  =  -jkO0  e 


J(©t-kx) 


and 


p  =  -jcop000  e 


J(cot-kx) 


(2.3) 

(2.4) 
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In  Equations  2.3  and  2.4,  <f>0  is  the  magnitude  of  the  velocity 
potential,  u  is  the  velocity  of  the  radiating  face,  to  is  the  angular 
frequency  and  k  is  the  wave  number. 

Recall  ZQ  is  the  acoustic  impedance  seen  by  a  transducer  when 

the  other  transducer  is  blocked.  Figure  3  represents  pictorially  this 
condition  (for  ease  of  analysis,  x=0  is  referenced  to  the  face  of  the 
blocked  driver  and  the  analysis  performed  for  the  left  driver).  The 
boundaiy  condition  at  x=0  requires  the  velocity  there  to  be  zero. 


Figure  3.  Blocked  Acoustic  Impedance. 

Applying  this  condition  results  in  the  following  forms  for  the  velocity 
potential,  acoustic  pressure  and  particle  speed  (the  periodic  time 


dependence  is  suppressed): 

O  =  O0  cos  (kx) .  (2.5) 

p  = -J(opo<I>o  cos  (kx) .  (2.6) 

u  = -kO0  sin  (kx).  (2.7) 

At  x=  -L  these  become: 

O  =  <t>0  cos  (kL) .  (2.8) 

p  = -jospo^ocos  (kL).  (2.9) 

u  =  k  O0  sin  (kL) .  (2.10) 
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For  the  transducer  at  x=-L, 

:o; 


Zn=^=  J- 


(2.11) 


U  uS 

where  S  is  the  cross  sectional  area  of  the  transducer  face. 
Substituting  Equations  2.9  and  2.10  into  Equation  2.11  gives 


~  _  -jco  Po *1*0  cos  (kL) 
0  S  k  <I>o  sin  (kL) 

Using  the  relationship  kc=to  yields 


(2.12) 


Z0  =  cot(kL).  (2.13) 

This  expression  for  ZQ  will  be  simplified  later  when  the  long 
wavelength  approximation  is  applied. 

3.  Transfer  Impedance 

T  is  a  transfer  impedance  and  is  the  ratio  of  the  pressure  on 
the  face  of  a  blocked  transducer  at  one  end  of  the  chamber  to  the 
volume  velocity  of  a  moving  transducer  at  the  other  end  (again,  x=0  is 
referenced  to  the  blocked  transducer  and  analysis  performed  on  the 
left  driver).  Figure  4  below  represents  this  situation. 


x - > 


S3  Blocked  1  I  Moving 


Figure  4.  Transfer  Impedance. 
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The  boundary  condition  at  x=0  requires  the  velocity  there  to  be 
zero.  Applying  this  condition  results  in  the  forms  expressed  in 
Equations  2.5  through  2.7.  Pressure  at  x=0  and  the  particle  velocity 
at  x=L  can  then  be  expressed  as 

p(0)  =  -j<op0<I>o  .  (2.14) 

and  u(L)  =  -k  4»0sin  (kL)  .  (2.15) 

By  inspection  of  Figure  2,  note  U2  =  -  S  u(L),  so  that 

fsijj°i=  Pffl  .  (2.16) 

U2  -S  u(L) 

Substituting  Equations  2.14  and  2.15  into  Equation  2.16  yields 

-jwpo^o 


T  = 


(2.17) 


S  k  <X>0  sin  (kL) 

This  can  further  be  simplified  to  the  following  expression: 

-JPoC 


T  = 


(2.18) 


S  sin  (kL) 

The  long  wavelength  approximation  will  now  be  applied  to  Z0  and  T. 

4.  Long  Wavelength  Approximation 

This  experiment  will  be  conducted  at  frequencies  for  which  the 
wavelength  is  much  larger  than  the  size  of  the  cavity.  The  sine  and 
cosine  terms  are  expanded  in  Taylor  series  form,  and  since  the 
product  (kL)  is  small,  the  terms  of  fourth  order  and  higher  are 
negligible.  The  sine  and  cosine  functions  are  expressed  as: 

,3 

(2.19) 


sin(kL)  =  kL 


cos(kL)  a  1 


(kLp 


3! 

2 


(kL) 

2! 


(2.20) 
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Applying  these  approximations  to  Equation  2.13  give 


^  _-JPoC  -J^Po 

2q  —  —  cot  (kL)  =  - 


S  k2L 


1  - 


(kL)- 


(kL|- 


Using  a  Taylor  expansion  on  the  denominator  i.e., 


-.-l 


1  - 


(kL)- 


=  1  + 


(kL)- 


the  above  equation  is  reduced  to  Equation  2.22: 

i2 


Z0  = 


-JPoc 


<kL)- 


kL 


(2.21) 


(2.21a) 


(2.22) 


Recalling  o>=kc  and  B  (bulk  modulus)  =  pc2. 


z0=-JL  +  «l. 

0  JcoSL  3S 


(2.23) 


Similar  relationships  are  applied  to  the  expression  for  T  found  in 
Equation  2.18,  resulting  in: 


f  =  B  Joy1' 

~  jcoSL  6S  ' 


(2.24) 


With  these  expressions,  the  acoustic  impedance  presented  to 
the  face  of  each  transducer  can  now  be  determined.  The  next  step  is 
to  apply  these  equations  to  the  case  when  the  two  transducers  are 
identically  driven  in-phase  and  out-of-phase. 
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5.  Identical  Transducers 

When  matched  transducers  are  used  in  the  cavity  and  are  driven 
identically  either  in-phase  or  out-of-phase,  a  solution  for  the  acoustic 
impedance  (Zac)  presented  to  each  transducer  face  is  readily 

obtained.  -Consider  Figure  2.  Using  this  figure,  the  acoustic  pressure 
at  the  face  for  each  transducer  can  be  expressed  as: 


Pi  ^ZqUj  +TU2.  (2.25) 

j^fUj  +Z0U2.  (2.26) 

If  the  transducers  are  identical  and  driven  in  phase,  U  i  =  U  2  =  U  . 

Equations  2.21  and  2.22  become: 

Pi  =  ZQU  +  TU  e  p.  (2.27) 

p2*fu  +  Z0Usp.  (2.28) 

Selecting  either  equation,  and  solving  for  the  ratio  of  acoustic  pressure 
to  volume  velocity  gives 

Zac  ~  Z0  +  T  ,  (2.29) 


where  Zac  is  the  acoustic  impedance.  Substituting  Equations  2.23  and 
2.24  into  the  above  gives 


Zor  = 


(2.30) 


'ac  "  jo>S(L/2)  ' 

This  is  the  same  as  for  a  closed  cavity  of  length  L/2.  In  this  case,  a 
pressure  antinode  exists  at  the  midplane  of  the  cavity. 

If  the  same  identical  drivers  are  driven  out  of  phase,  the  same 
approach  is  taken  and  the  result  is: 


^  _jcop(L/2) 
^ac=  c 


(2.31) 
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This  is  the  same  as  for  an  open  cavity  of  length  L/2  (neglecting 
radiation).  In  this  case,  a  pressure  node  exists  at  the  midplane  of  the 
cavity. 


6.  Fluid-Filled  Porous  Solid 

Equations  have  been  derived  for  the  acoustic  impedance  of  a 
chamber  filled  only  with  fluid.  When  sound  propagates  in  a  fluid  - 
filled,  rigid  porous  body,  acoustic  energy  losses  occur  due  to  the 
thermal  conductivity  and  shear  viscosity  of  the  fluid.  Consequently, 
the  bulk  modulus  and  effective  mass  density  become  complex. 
Equations  2.30  and  2.31  now  become: 


B 

j<nPS(L/2) 


,  and 


(2.32) 


Z 


_  j©(L/2)  ~ 
ac~  ps  p  • 


(2.33) 


The  additional  parameter,  P,  is  the  porosity.  The  porosity  is  defined 
as  the  ratio  of  the  volume  occupied  by  the  fluid  to  the  total  sample 
volume.  It  is  required  since  the  cross  sectional  area  for  fluid 
interaction  is  reduced  by  the  frame  of  the  porous  material.  A  formal 
derivation  of  the  porous  media  effects  is  presented  in  the  Appendix. 

C.  TRANSDUCERS 

Since  this  research  involves  determining  complex  mechanical 
quantities  from  electrical  impedance  measurements,  it  is  important  to 
first  identify  the  impedance  of  the  transducers  used.  The  transducers 
used  are  antireciprocal,  moving  coil  type.  The  theory  for  this  type  of 
transducer  is  well  understood  (Ref.  81. 


Equations  have  been  developed  to  explain  the  relationships 
between  voltage  applied  to  a  transducer  and  the  force  on  its  radiating 
face.  The  equations  are: 

V  =  ZEB  I  +  (Bl)  u,  (2.34) 

F  =  -(Bl)  I  +  ZnjQ  u.  (2.35) 

In  the  above  equations: 

V  =  the  voltage  across  the  input  of  the  transducer. 

I  =  the  input  current  to  the  transducer. 

F  =  the  force  on  the  radiating  surface, 
u  =  the  velocity  of  the  radiating  surface. 

ZEB  =  the  blocked  electrical  impedance  (V/f;  u=0). 

Zjjjq  =  the  open  circuit  mechanical  impedance  (F/u;f=0). 

(Bl)  =  the  product  of  the  magnetic  field  and  length  of  coil 
in  the  gap  of  the  driver  =  Transduction  coefficient. 

The  moving  coil  transducer  can  further  be  modeled  by  a  simple 

equivalent  electrical  circuit  (in  the  mobility  analog)  given  in  Figure  5 

(Ref.  7:p.  352): 

Ro  L0 

/V 

V 

Figure  5.  Driver  Equivalent  Circuit. 

In  the  figure  above, 

Ro  =  Blocked  electrical  resistance, 

Lq  =  Blocked  electrical  inductance, 

Rm  =  Electrical  equivalent  mechanical  resistance, 
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Lm  =  Electrical  equivalent  mechanical  inductance, 

CM  =  Electrical  equivalent  mechanical  capacitance. 

The  electrical  equivalents  are 

obtained  from  the  following 

relationships  [Ref.  9): 

R  -  (B1)2 

RM  -  R 

I'-mech 

(2.36) 

L  -  (B1)2 

Lm-  k  • 

(2.37) 

_  m 

Cm-  2  • 

(2.38) 

(Bl) 

The  parameters  RmeCh*  k  and  m  are  the  mechanical  resistance, 
stiffness  and  mass  of  the  driver,  respectively.  From  the  above 
relationships,  if  the  stiffness,  moving  mass,  mechanical  resistance, 
transduction  coefficient  and  blocked  electrical  impedance  are  known, 
the  electrical  impedance  of  the  driver  as  a  function  of  frequency  can 
be  predicted. 

Additionally,  for  the  transducers  selected,  there  is  an  extra 
stiffness  existing  due  to  the  air  in  the  cavity  behind  the  dome  of  the 
voice  coil.  This  stiffness  should  exhibit  the  characteristics  of  an 
acoustic  cavity  [Ref.  10)  given  by, 

Y  S2P0 

k  = - - - ,  (2.39) 

where 

Y  =  ratio  of  the  constant  pressure  to  constant  volume 
specific  heat  capacity  for  a  gas, 
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S  =  the  area  of  the  driving  face, 

PG  =  the  ambient  pressure, 

V  =  the  cavity  volume. 

This  added  stiffness  modifies  the  equivalent  electrical  circuit 
presented  in  Figure  5.  The  modified  circuit  is  illustrated  in  Figure  6. 
The  new  inductance,  ,  accounts  for  the  additional  dynamic 

stiffness.  Note  this  model  does  not  reflect  any  acoustic  loading  on  the 
transducer.  The  effects  of  acoustic  loading  on  the  electrical 
properties  of  the  transducers  will  be  presented  in  the  next  section  of 
this  chapter. 


Figure  6.  Modified  Driver  Equivalent  Circuit. 


D.  COUPLER-DRIVER  SYSTEM 

When  the  electrical  and  mechanical  characteristics  of  the  two 
drivers  are  combined  with  the  acoustical  characteristics  of  the 
chamber,  the  two  port  system  is  given  in  Figure  7.  In  the  figure 
below,  the  Zr  are  the  electrical  equivalent  impedances  of  the  radiation 

impedances  presented  to  the  transducers.  The  radiation  impedance 
as  used  here  is  a  mechanical  impedance,  i.e.,  a  ratio  of  force  to 
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Ze,  2e2 


Figure  7.  Combined  Electrical  Circuit. 


velocity.  This  differs  from  the  acoustic  Impedance  by  a  factor  of  S2, 
where  S  Is  the  moving  area  of  the  driver. 

When  the  drivers  are  driven  identically  in  parallel  "in-phase"  or 
"out-of-phase".  Figure  7  reduces  to  the  equivalent  circuit  in  Figure  8, 
where  I  /2  is  the  current  to  flow  through  each  driver  and  V  is  their 
common  voltage. 


Ze 


T 

V 


Figure  8.  Combined  Electrical  Circuit  For  Identical  Drivers. 


Therefore,  by  measuring  the  input  electrical  impedance  of  the  parallel 
combination  (V/I),  and  knowing  Ze  and  ZM,  Zr  can  be  determined, 
and  from  it,  the  value  of  Zac  extracted  according  to: 


zin= 


V 

I 


1 

2 


Ze  + 


ZmZr 
Zm  +  Zr  , 


(2.40) 
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Defining  Z  ■  2  Z 


IN 


Zc-  = 


ZM  ZR 
ZM  +  ZR 


then 


Zp  = 


ZMZ 
ZM  "  Z 


(2.41) 

(2.42) 


The  acoustic  impedance  is  extracted  from  ZR  in  the  mobility  analog 


by. 


zac  = 


2  2 
(Bl)  S 

ZR 


(2.43) 


Once  Zac  is  determined,  the  complex  bulk  modulus  is  extracted 
using  Equation  2.32  when  the  drivers  are  in-phase: 

B  =  Jffl  (t)  PS  Zac  ,  (2.44) 

or  the  complex  mass  density  using  Equation  2.33  when  out-of- phase: 

PS  ^ 


P= 


(2.45) 


J<o  (L/2)  ac' 

In  summary,  by  measuring  the  input  electrical  impedance  of  two 
identical  drivers  used  as  the  end  caps  of  a  small  cylindrical  chamber, 
the  complex  parameters  B  and  p  of  the  fluid  in  a  porous  rigid  material 
can  be  determined,  provided  the  characteristics  of  the  drivers  are 
known. 
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HI.  TRANSDUCER  CHARACTERIZATION  AND  COUPLER  DESIGN 

A.  BACKGROUND 

The  experimental  portion  of  this  investigation  was  conducted  in 
two  major  phases.  The  first  phase  involved  the  selection  and 
characterization  of  two  matched  transducers.  Since  it  was  crucial  to 
characterize  the  transducers  as  accurately  as  possible,  several 
experiments  were  conducted,  and  where  possible,  alternative 
methods  used  to  verify  previous  results. 

The  second  major  phase  involved  the  design  and  construction  of  a 
test  chamber  (coupler).  The  coupler  was  used  to  hold  the  selected 
transducers  and  the  test  sample.  Once  this  was  completed, 
measurements  were  made  to  verify  the  use  of  the  input  electrical 
impedance  to  extract  physical  parameters. 

In  this  chapter  the  experimental  apparatus  used  in  each  phase  is 
described.  Accompanying  the  description  is  an  explanation  of  how  the 
apparatus  functions.  Finally,  sample  data  are  presented. 

B.  TRANSDUCER  SELECTION  AND  PARAMETER 

CHARACTERIZATION 

1.  Matched  Transducer  Selection 

Matched  transducers  can  be  obtained  by  either  constructing  the 
transducers  from  identical  components  or  by  obtaining  many  samples 
of  a  manufactured  product  and  testing  several  units  until  two  matched 
transducers  are  found.  The  latter  method  was  selected,  primarily  due 
to  the  experience  of  Hofler  [Ref.  11].  Hofler  had  worked  with 
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acoustically  "soft"  (moving  coil)  transducers.  The  type  of  transducer 
he  used  would  meet  the  requirements  of  this  experiment. 

The  transducers  are  nominal  one-inch  diameter  domed 
tweeters.  They  are  manufactured  by  Philips,  Model  AD062T8 
[Ref.  12).  The  dome  is  made  of  a  water-resistant  polyvinyl  material. 
The  cover  of  the  driver  is  easily  removed.  Once  this  is  done,  the  voice 
coil  and  its  support  ring  can  be  readily  detached  from  the  magnet.  A 
cutaway  view  of  the  driver  (minus  cover)  is  provided  in  Figure  9. 


Voice  Coil 

Figure  9.  Cutaway  View  of  Philips  Driver. 


Several  drivers  were  obtained  and  the  impedance  as  a  function 
of  frequency  determined  using  a  Hewlett-Packard  Impedance 
Analyzer,  Model  4194A  (HP-4194A)  [Ref.  13],  as  follows.  The  drivers 
were  mounted  on  a  large  baffle  and  placed  in  an  anechoic  chamber. 
The  resistive  and  reactive  components  of  the  impedance  of  each 
driver  were  measured  over  the  frequency  range  from  100  Hz  to  10 
kHz.  The  impedance  curves  for  each  driver  were  then  compared  and 
two  drivers  with  similar  responses  were  selected  for  further  testing. 
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Table  3.1  presents  the  resistance  and  reactance  values  of  the 
two  transducers  selected  at  their  respective  mechanical  resonance 
frequencies.  The  similarities  in  resonance  values  resulted  in  their 
selection  for  further  experimentation. 


|  TABLE  3.1.  RESONANCE  VALUES  FOR  SELECT! 

ED  DRIVERS. 

Driver 

Frequency  (Hz) 

Resistance(ohms) 

Reactance  (ohms) 

1 

954.0 

65.19 

1.36 

2 

954.0 

70.31 

4.36 

2.  Source  Strength 

To  further  confirm  the  similarity  of  the  two  drivers,  the  source 
strength  of  each  was  found.  The  source  strength  measurement  is  best 
explained  by  the  use  of  Figure  10. 

The  HP-4194A  provides  the  current  for  the  driver.  In  the  Gain- 
Phase  mode,  it  is  able  to  measure  the  voltage  difference  between  a  test 
and  a  reference  channel,  or  display  the  voltage  of  each  channel.  In 
this  case,  the  reference  is  the  voltage  associated  with  the  current 
through  a  10  ohm  resistor  on  the  return  side  of  the  driver.  The 
HP-4194A  was  swept  through  the  desired  frequency  span.  The  output 
was  then  printed  on  an  accompanying  printer  and  stored  in  the 
impedance  analyzer  for  further  analysis.  The  distance  was  measured 
from  speaker  to  microphone  for  source  strength  calculations. 

The  calibration  standard  was  a  General  Radio  (GR)  microphone. 
The  calibration  of  this  microphone  was  accomplished  by  use  of  a 
pistonphone.  The  calibration  was  further  checked  by  using  a  less 
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accurate  GR  Calibrator.  The  sensitivity  level  obtained  by  these  two 
methods  was  then  verified  against  the  calibration  curves  accompanying 
the  microphone. 


Sweeps  were  conducted  for  the  two  selected  drivers.  Printouts 
of  the  reference  and  test  channels  were  obtained  for  each  driver,  and 
then  compared.  When  the  output  from  the  two  different  drivers  were 
compared,  they  were  virtually  identical.  Based  upon  this,  these 
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transducers  were  retained  for  the  remaining  portions  of  the 
experiment. 

3.  Transducer  Parameter  Determination 

Five  parameters  must  be  known  to  adequately  predict  the 
performance  of  the  moving  coil  transducers.  These  parameters  are 
the  blocked  electrical  impedance  (ZBE).  driver  stiffness  (k),  moving 

mass  (m),  the  transduction  coefficient  (Bl)  and  mechanical  resistance 
(Rm)* 

a.  Electrical  Parameters 

The  easiest  electrical  parameters  to  determine  are  the 
blocked  electrical  parameters.  Two  methods  were  employed  to 
accomplish  this.  In  the  first  method  the  voice  coil  was  removed  from 
the  driver.  With  the  coil  removed,  a  digital  multimeter  was  used  to 
determine  the  direct  current  resistance  of  the  coil.  By  use  of  a 
Hewlett-Packard  Low  Frequency  Impedance  Analyzer,  Model  4192 
(HP-4192),  the  inductance  at  several  frequencies  was  measured  and 
then  averaged.  The  second  method  was  similar  to  the  first  except  the 
voice  coil  remained  in  the  driver.  The  moving  portion  was  held  fast. 
Resistance  and  inductance  were  measured  over  several  frequencies 
and  averaged. 

The  impedance  obtained  by  these  two  methods  was  used  to 
determine  the  best  value  for  comparing  modeled  to  measured 
performance.  Table  3.2  presents  typical  data  obtained  by  these 
measurements. 
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The  DC  resistance  did  not  vary  appreciable  from  the  in¬ 
magnet  case  to  the  out-of-magnet  case.  However,  the  blocked 
inductance  changed  by  a  factor  of  ten  from  the  in-magnet  to  out-of¬ 
magnet  case.  The  out-of-magnet  value  was  selected  for  use  since  it 
provided  the  best  fit  to  the  modeled  performance  of  the  lumped 
parameters.  The  in-magnet  inductances  were  difficult  to  obtain  since 
it  was  a  dynamic  measurement  and  blocking  the  voice  coil  was  not 
easy. 


TABLE  3.2.  BLOCKED  ELECTRICAL  PARAMETERS. 

|  Blocked  Resistance  (ohms) 

Blocked  Inductance  (mH)  | 

Coil  Removed 
Driver 

1 

6.544 

0.11  (at  30  Hz) 

2 

6.542 

0.11  (at  30  Hz) 

Coil  Installed 
Driver 

1 

6.546 

1.07  (at  30  Hz) 

2 

6.508 

1.10  (at  30  Hz) 

Since  it  was  crucial  to  know  the  blocked  impedance  as 
accurately  as  possible,  the  diaphragm  of  a  similar  driver  with  identical 
out-of-magnet  electrical  parameters  was  rigidly  clamped  to  the 
magnetic  structure  using  epoxy  resin.  When  the  epoxy  was  cured,  the 
blocked  inductance  was  measured  as  a  function  of  frequency.  The 
results  are  presented  in  Figure  11.  The  data  reveals  the  blocked 
inductance  is  frequency  dependent.  To  properly  account  for  this 
dependence,  the  blocked  impedance  over  the  desired  frequency  range 
will  be  obtained  from  the  epoxied  driver.  The  frequency  dependent 
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Blocked  Inductance  (jiHenries) 


blocked  impedance  data  can  then  be  stored  in  the  data  registers  of  the 
HP-4194A  for  later  use  when  the  acoustic  impedance  experiments  are 
conducted. 


Blocked  Inductance 


b.  Driver  Stiffness 

Several  methods  can  be  used  to  determine  the  stiffness  of 
the  transducers.  The  first  method  attempted  was  a  static  method. 

i 
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The  motivation  for  making  a  static  determination  is  the  ease  of 
measurement.  If  a  known  mass  Am  is  added  to  the  dome  and  the 
voice  coil  deflection  Ax  is  measured,  the  stiffness  can  be  determined 
by: 

Amg  =  kAx.  (3.1) 

In  the  experiment  performed,  the  voice  coil  deflection  was 
measured  with  a  linear  variable  differential  transformer  (LVDT).  The 
LVDT  operates  with  alternating  current  (ac)  and  has  the  property  that 
the  voltage  across  the  secondary  is  linearly  proportional  to  the 
position  of  the  plunger.  If  the  change  in  LVDT  secondary  voltage  is 
measured  and  the  LVDT  calibrated,  Equation  3. 1  can  be  written  as: 

Amg=kAV~  .  (3.2) 

where  AV  is  the  output  voltage  of  the  LVDT  and  dx/dV  is  the 
reciprocal  of  the  slope  of  the  curve  (line)  of  output  voltage  versus 
displacement. 

The  LVDT  used  was  manufactured  by  the  Collins  Corporation, 
Model  LMS-1 19V22-01  [Ref.  14J.  It  was  calibrated  using  a  calibration 
stand  and  a  micrometer.  The  calibration  stand  held  the  transformer 
securely  in  place  while  the  transformer's  lightweight  plunger  was 
allowed  to  rest  atop  the  micrometer.  As  the  micrometer  position  was 
changed,  the  output  of  the  LVDT  secondary  was  recorded.  The  data 
obtained  were  plotted  as  LVDT  secondary  output  voltage  versus 
micrometer  position  (Figure  12).  The  slope  of  this  line  represents 
the  change  in  voltage  for  a  change  in  position.  Table  3.3  lists  the 
slope  and  its  relative  uncertainty. 
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LVDT  Calibration  Curve 


Position  (in) 

Figure  12.  LVDT  Calibration  Curve. 


TABLE  3.3.  LVD1 

r  CALIBRATION. 

Slope  [mV  (rms)  /  in] 

Relative  Uncertainty  {%) 

|  -996.2 

0.19 

A  special  ring  was  machined  from  lightweight  aluminum  to 
fit  smoothly  around  the  voice  coil  dome.  Ring  shaped  brass  weights 
were  added  to  the  special  ring.  A  1000  Hz  stabilized  sine  wave  was 
applied  to  the  primaiy  side  of  the  LVDT,  whose  plunger  rested  upon 
the  top  of  the  voice  coil  dome.  The  change  in  LVDT  secondary  voltage 
caused  by  the  change  in  position  due  to  the  weight  of  the  brass  ring 
was  measured  using  a  Stanford  Research  Lock-in  Analyzer,  Model  530 
(Ref.  15).  Table  3.4  lists  the  data  points  obtained  using  this  method. 
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TABLE  3.4.  STATIC  STIFFNESS  DATA. 


1 -  1 

|  Mass  Added  (g) 

LVDT  Output  (mV)  | 

|  Driver  1  14.98 

2.663 

41.80 

6.751 

84.84 

13.26 

109.90 

17.45 

Driver  2  14.98 

2.693 

40.49 

6.714 

84.84 

14.52 

109.90 

17.35 

Two  methods  can  be  used  to  reduce  the  stiffness  data.  Both 
were  utilized  to  provide  a  cross-check.  The  first  involves  plotting  the 
change  in  LVDT  secondary  voltage  as  a  function  of  mass  added  to  the 
driver.  The  slope  of  this  line  is: 

slope=1kS-  <3-3' 


Solving  for  k: 

k  = 


g  dV 


(3.4) 


Slope  dx 

The  second  method  of  reduction  involves  calculating  k 
directly  from  each  data  point: 


AmgdV 

AV  dx  ‘  l3’5) 

These  values  for  stiffness  are  then  plotted  versus  the  force  applied. 
The  result  should  be  a  nearly  horizontal  line  whose  intercept  reflects 
the  value  of  stiffness  for  small  displacements. 

Both  schemes  were  employed  to  analyze  the  data.  Due  to  a 
slight  but  observable  nonlinearity  in  the  results  from  the  first  method, 
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the  second  method  was  selected  for  final  data  reduction.  The 
transducers  are  less  stiff  at  low  displacements  and  tend  to  be  more 
stiff  at  high  displacements.  Since  operation  of  the  driver  in  the 
chamber  experiment  will  be  in  the  low  displacement  region,  it  was 
believed  method  two  would  provide  a  more  realistic  value  for  stiffness. 

Table  3.5  lists  the  stiffness  and  relative  uncertainty  for  each 
method  of  reduction.  The  difference  due  to  the  slight  nonlinearity  is 
seen  by  the  difference  in  values. 


TABLE  3.5  STIFFNESS  OF  SELECTED  DRIVERS. 

Reduction  Method 

Stiffness  (N/m) 

Relative  Uncertainty  (%) 

One 

Driver 

1  1 

2448 

1.6 

2 

2380 

2.9 

Two 

Driver 

1 

2225 

4.3 

2 

2118 

4.4 

Another  way  to  determine  stiffness  is  to  extract  it  from  the 
resonance  frequency  (co=2xf)  and  the  (presumed  known)  moving  mass: 


Since  to  and  m  are  known,  this  equation  can  be  rearranged  to  solve  for 
the  stiffness,  k: 

k  =  m(2xf)2  .  (3.7) 

This  was  the  method  finally  selected  because  it  eliminated  a  sizeable 
discrepancy  that  existed  when  the  modeled  electrical  impedance  of  a 
driver  was  compared  with  its  measured  electrical  impedance.  As  will 
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be  discussed  in  the  next  section,  the  moving  mass  of  an  auxiliary 
driver  was  measured.  Using  this  mass  and  the  vacuum  resonance 
frequency  of  each  of  the  selected  drivers,  the  stiffness  for  each 
selected  driver  was  solved  for  and  is  listed  in  Table  3.6.  The  relative 
uncertainty  is  derived  from  the  uncertainties  of  the  moving  mass  and 
the  resonance  frequencies. 


I-  TABLE  3.6.  STIFFNESS  OF  SELECTED  DRIVERS  I 

Driver 

Stiffness  (N/m) 

Relative  Uncertainty  (%) 

1 

5089 

3.34 

2 

4687 

3.34 

c.  Moving  Mass 

Several  methods  can  be  used  to  determine  the  moving  mass 
of  the  transducers.  The  first  is  the  method  of  added  mass.  With  the 
driver  operating  in  a  vacuum,  the  resonance  frequency  is  given  by: 


where  k  is  the  stiffness  constant  and  m  the  moving  mass  of  the  driver. 
If  a  mass  mj  is  added  to  the  voice  coil  dome,  the  resonance  frequency 

now  becomes: 

■  (3-9) 

By  rearranging  this  formula,  plotting  the  square  of  the  period  as  a 
function  of  mass  added,  the  slope  of  the  line  is  4 n2  /  k  and  the 
intercept  is  (4rt2  m)  /  k.  Dividing  the  intercept  by  the  slope  will  yield 
the  moving  mass  of  the  driver. 


28 


The  method  of  using  added  masses  to  determine  the  moving 
mass  of  the  driver  was  abandoned  because  of  inconsistent  results. 
This  was  due  apparently  in  the  inability  to  obtain  satisfactorily  stiff 
bonding  of  the  added  mass  at  frequencies  of  interest. 

A  second  method  is  very  simple.  Assuming  all  the  drivers 
are  identical  in  construction,  the  mass  of  any  given  driver  should  equal 
that  of  another,  within  the  tolerance  of  manufacturing.  Using  this 
assumption,  an  estimate  of  the  moving  mass  can  be  obtained  by 
measuring  the  mass  of  a  diaphragm  removed  from  an  extra  driver. 

This  was  accomplished  as  follows.  The  diaphragm  of  an 
extra  driver  was  cut  away  from  the  support  ring.  The  surround  was 
then  removed  from  the  diaphragm.  The  mass  of  each  was  measured 
separately,  the  results  shown  in  Table  3.7. 


s  TABLE  3.7  MOVING  MASS  OF  SELECTED  DRIVERS.  | 

Mass  of  Dome  and  Voice  Coil(g) 

Mass  of  Surround  (g) 

0.3086 

0.0286 

An  initial  estimate  of  the  moving  mass  was  taken  as  the  mass 
of  the  voice  coil  and  its  diaphragm  plus  one-half  the  mass  of  the 
surround.  The  value  was  0.2943  grams,  with  a  relative  uncertainty  of 
3  per  cent  based  on  half  the  mass  of  the  surround.  This  value  was 
later  refined  by  curve  fitting  input  electrical  impedance  data  to  an 
equivalent  circuit  model,  described  in  Section  II. C. 
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d.  Transduction  Coefficient 


The  transduction  coefficient  (Bl)  is  obtained  directly  as 
described  below.  The  experimental  setup  is  the  same  as  the  stiffness 
experiment.  The  two  experiments  are  compatible  and  can  be 
conducted  together. 

The  driver  was  placed  on  a  table  facing  up.  A  mass  was 
added  to  the  dome  of  the  voice  coil,  resulting  in  a  certain 
displacement  of  the  diaphragm.  A  direct  current  was  applied  to  the 
driver  and  adjusted  to  return  the  diaphragm  to  its  original  rest 
position,  as  indicated  by  the  output  voltage  of  a  LVDT.  By  knowing  the 
mass  added,  the  applied  force  was  calculated  from  F=mg.  For  this 
known  force,  the  current  required  to  exactly  offset  it  was  measured. 
Using  the  relationship  F=mg=Bli,  the  Bl  product  was  found.  Several 
masses  were  added  and  the  Bl  product  calculated  for  each.  The  Bl 
products  were  then  averaged  and  the  standard  deviation  calculated. 
Table  3.8  shows  the  data  for  Driver  1  obtained  this  way. 


TABLE  3.8  TRANSDUCTION  COEFFICIENT  DATA. 

Mass  Added  (g) 

Current  (ma  DC) 

Transduction  Coef.  (T-m) 

14.98 

44.93 

3.268 

28.40 

84.46 

3.295 

41.80 

126.2 

3.246 

54.97 

173.5 

3.105 

67.30 

205.7 

3.207 

84.84 

256.9 

3.237 

l  l 

|  Average  Transduction  Coef.  (T-m) 

Relative  Uncertainty  (%)  1 

1  3.226 

2.0  I 
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A  second  method  was  used  to  check  the  values  for  the  B1 


product.  Using  the  equation  F  =  Bli  =  -kx,  the  derivative  of  this 
equation  is  taken  with  respect  to  current.  This  results  in  the 
expression: 

Bl=-k^.  (3.10) 

dl 

The  stiffness  (k)  is  obtained  from  by  the  procedure  listed  in  Section 
III.B.3.b. 

The  term  (dx/dl)  is  determined  by  applying  a  known  direct 
current  to  the  speaker  and  measuring  the  LVDT  output.  Table  3.9 
lists  sample  data  for  Driver  1. 


TABLE  3.9.  CURRENT  RESPONSE. 

Current  (ma  DC) 

LVDT  Output  (mV) 

-251.02 

307.86 

148.81 

302.57 

-93.79 

299.47 

56.26 

297.19 

-2.92 

294.08 

46.42 

291.79 

j  102.82 

288.30 

144.66 

285.10 

256.11 

276.60 

Linear  regression  was  performed  on  the  above  data,  with  the  current 
the  independent  variable,  and  the  LVDT  output  voltage  the  dependent 
variable.  The  slope  obtained  by  this  linear  regression  is  dV/dl.  For 
Driver  1,  the  change  in  LVDT  output  with  respect  to  direct  current 
applied  was  0.05609  mV/ma.  Knowing  this  and  the  LVDT  calibration. 


Applying  Equation  3.5  to  Equation  3.11,  the  result  is: 


Am  g  dV  .  (3.12) 

(AV)  dl 

The  transduction  coefficient  was  calculated  for  each  mass 
added,  using  the  respective  change  in  LVDT  output.  The  output  of  the 
LVDT  as  a  function  of  current  applied  was  determined  as  discussed 
above.  The  B1  values  were  averaged  and  the  standard  deviation  was 
used  to  calculate  relative  error.  This  resulted  in  the  value  of  B1  listed 
in  Table  3. 10.  The  value  is  slightly  higher  than  the  one  listed  in  Table 
3.8.  This  was  a  result  of  the  slight  nonlinearity  in  stiffness.  The  first 
method  used  essentially  a  lower  value  of  stiffness  (displacement  was 
always  near  zero).  The  second  method  incorporated  a  higher  value  of 
stiffness  since  the  stiffness  used  was  the  slope  of  the  force  versus 
displacement  for  a  wide  variety  of  displacements.  This  slope  was 
influenced  by  the  increasing  stiffness  with  increasing  displacement, 
hence  a  higher  transduction  coefficient  (displacement  was  over  the 
range  of  driver  motion). 


As  a  further  check  of  Bl,  a  measurement  of  the  B  field  in  the 
gap  was  made  with  a  Hall  probe.  The  number  of  turns  were  counted 
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with  a  microscope  and  the  length  of  wire  estimated  from  this  and  the 
diameter  of  the  voice  coil.  The  result  was  a  B1  product  of  2.8  T-m. 
This  aided  in  confirming  the  value  of  the  transduction  coefficient. 

e.  Mechanical  Resistance 

Mechanical  resistance  is  obtained  from  measured  values  of 
electrical  impedance  near  mechanical  resonance.  At  resonance  in 
vacuum,  the  magnitude  of  the  impedance  is  approximately  Rq  + 

(the  reactive  component,  jcoLo  ,  was  negligible)  from  the  equivalent 

circuit,  Figure  5.  Using  this  fact,  the  mechanical  resistance  is 

obtained  by  rewriting  Equation  2.36  in  the  form: 

R  -  (Bl)2  ,, 

Rm"  ~  |  p  (3.13) 

■^RES1 '  kO 


Table  3.11  lists  the  resonance  frequency,  resonance  impedance  and 
resulting  Rm  for  each  of  the  drivers. 


TABLE3.il.  VACUUM  RESONANCE  VALUES.  I 

Driver 

Res.  Freq.  (Hz) 

Impedance  (ohms) 

Mech  Resistance  (N-s/m) 

1 

662,00 

63.11 

0.1883 

2 

635.00 

50.71 

0.2172 

f.  Cavity  Stiffness 

From  the  cross  section  (Figure  9)  view  of  the  driver,  one  can 
observe  a  large  cavity  behind  the  dome  of  the  speaker.  This  cavity 
provides  significant  additional  stiffness  to  the  system  when  operated 
in  air.  Since  the  drivers  will  be  operated  at  various  static  pressures, 
the  stiffness  as  a  function  of  static  pressure  must  be  either  known  or 
eliminated. 


An  attempt  to  eliminate  the  stiffness  of  this  cavity  was  tried 
first.  A  small  hole  was  drilled  through  the  pole  piece  of  an  auxiliary 
driver.  The  input  electrical  characteristics  were  measured  for  several 
different  hole  diameters.  Instead  of  eliminating  the  stiffness,  the 
driver  now  had  additional  loading  caused  by  the  air  in  the  hole.  This 
resulted  in  a  coupled  oscillator  response  which  would  only  increase 
the  complexity  of  the  problem  of  added  stiffness. 

Since  the  cavity  stiffness  could  not  be  eliminated,  the  next 
step  was  to  determine  the  stiffness  as  a  function  of  pressure. 
Assuming  the  cavity  stiffness  is  due  only  to  the  gas  in  the  cavity,  it 
should  be  proportional  to  static  pressure  (see  Equation  2.39).  To 
determine  the  pressure  dependence  of  the  cavity  stiffness,  the  driver 
was  placed  in  a  pressure  chamber  and  the  pressure  was  changed  in 
controlled  increments.  Resonance  frequency  was  recorded  as  a 
function  of  pressure.  The  resonance  frequency  dependence  on  the 
cavity  stiffness  is  shown  in  Equation  3.14: 

f2  = - -  (3. 14) 

(4n  )(m  +  mr  ) 

where: 

f  =  measured  resonance  frequency. 

k  =  static  stiffness  of  the  driver. 

k(p)  =  the  pressure  dependent  cavity  stiffness. 

m  =  moving  mass  of  the  driver. 

mr  =  radiation  mass. 

Table  3.12  lists  the  data  obtained  from  this  experiment  using  Driver  1. 
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TABLE  3.12.  CAVITY  STIFFNESS  DATA.  H 

_ 

Pressure  (in  Hg) 

Res.  Frequency  (Hz) 

28.6 

650.00 

23.3 

727.99 

20.2 

772.00 

14.8 

854.00 

12.3 

888.00 

9.9 

912.00 

-0.51 

1016.99 

The  square  of  the  frequency  was  plotted  as  a  function  of 
pressure.  The  result  was  a  straight  line  (the  radiation  mass,  mr,  was 

negligible  for  the  drivers  at  these  frequencies).  Once  a  linear 
relationship  had  been  verified,  the  stiffness  as  a  function  pressure  was 
solved  for  at  each  frequency.  This  stiffness  was  then  plotted  as  a 
function  of  pressure.  The  result  was  another  straight  line,  expressing 
cavity  stiffness  as  a  function  of  pressure.  Figure  13  is  a  plot  of  these 
data.  The  equation  of  the  line  was  solved  for  by  regression,  resulting 
in: 

k(p)  =  (-159.4 +0.07350  (Pq))^-  .  (3.15) 

The  relative  error  for  a  value  from  this  equation  is  3.9%.  The 
intercept  value  for  this  equation  was  not  forced  to  be  zero.  The 
negative  value  is  a  result  of  the  random  errors  associated  with  this 
portion  of  the  experiment. 

To  check  the  validity  of  Equation  3.15,  the  entire  volume 
behind  the  dome  of  an  auxiliary  speaker  was  measured  by  filling  it  with 
a  known  volume  of  fluid.  Using  Equation  2.39,  the  stiffness  was 
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calculated  for  a  pressure  equal  to  one  atmosphere.  The  value  agreed 
to  within  one  percent  of  the  value  given  by  Equation  3. 15. 

g.  Modeled  Performance 

Once  all  the  electrical  and  mechanical  parameters  were 
identified  for  each  driver,  the  parameters  were  converted  to  their 
equivalent  electrical  impedances  using  Equations  2.36  through  2.38. 
Utilizing  a  circuit  analysis  program,  the  impedance  of  the  circuit  in 
Figure  5  was  modeled  on  a  computer.  The  resulting  modeled 
electrical  impedance  values  were  then  compared  to  the  actual 

Cavity  Stiffness  as  a  Function  of  Pressure 


Figure  13.  Cavity  Stiffness  Pressure  Dependence, 
electrical  impedance  measurements  made  earlier.  The  first  model 
incorporated  the  electrical  equivalent  for  stiffness  obtained  using 
added  mass  (Section  III.B.3.b),  and  the  electrical  equivalent  for  mass 
extracted  from  resonance  frequency  (Equation  3.8).  The  results  of 
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this  model  are  given  in  Figure  14  (only  the  magnitude  of  the 
impedance  is  shown). 

To  resolve  this  discrepancy  in  the  fit,  the  moving  mass  of  an 
auxiliary  transducer  was  measured  and  stiffness  obtained  from  the 
resonance  frequency  (as  described  in  Sections  III.B.3.b  and  III.B.3.C). 
Figure  15  demonstrates  the  performance  based  upon  the  revised 
model.  The  fit  has  improved  significantly. 

70 

60 

1 

O  50 

C  40 
m 
■a 

f  » 

20 

10 

0 

100  1000  10000 

Frequency  (Hz) 

Figure  14.  Modeled  Performance  With  Original  Mass  Estimate. 

A  best  fit  was  attempted  using  the  built-in  circuit  analysis 
program  of  the  HP-4194A.  This  further  refined  the  values  used  in  the 
model  which  will  be  used  to  extract  the  acoustic  impedance  from  the 
input  electrical  impedance.  Figure  16  demonstrates  the  excellent  fit 


Modeled  vs  Measured  using  Original  Data 
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of  measured  versus  modeled.  A  best  fit  was  performed  for  each  driver. 
The  values  were  then  averaged  to  define  a  single  transducer  whose 
characteristics  would  be  used  for  the  "identical  transducfers".  Table 
3.13  fists  the  specific  parameters  and  their  electrical  equivalents  to 
be  used  in  the  circuit  model  (Figure  6)  to  extract  the  value  of  acoustic 
impedance. 

Modeled  vs  Measured  Using  Removed  Mass 


100  1000  10000 
Frequency  (Hz) 


Figure  15.  Modeled  Performance  With  Removed  Mass  Estimate. 
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TABLE  3.13.  MODEL  TRANSDUCER  PARAMETERS 

II  Parameter 

Value 

Units  || 

|  Blocked  Resistance 

6.3521 

ohms  H 

|  Blocked  Inductance 

0.11 

mH  | 

|  Mechanical  Resistance 

49.08 

ohms  1 

|  Driver  Mass 

26.65 

HF 

1  Driver  Stiffness 

2.269 

mH  1 

H  Cavity  Stiffness 

1.351 

mH  1 

Modeled  vs  Measured  for  Revised  Circuit  Elements 
70 

60 

I 

o  50 

40 

30 

20 

10 

0 

100  1000  10000 

Frequency  (Hz) 

Figure  16.  Best  Fit  Model. 
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C.  COUPLER  DESIGN  AND  CONSTRUCTION 
1.  Coupler  Construction 

Figure  17  is  a  cutaway  view  of  the  assembled  cylindrical  brass 
coupler  (for  clarity,  only  some  dimensions  are  shown).  The  main 
parts  of  the  coupler  are  the  center  chamber  with  sample  sleeve  and 
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two  hermetically  sealed  driver  chambers,  one  on  each  end.  The 
drivers  are  secured  in  their  respective  chambers  with  a  backing  plate. 

The  dimensions  of  the  coupler  were  dictated  in  part  by  the 
dimensions  of  the  driver  and  by  the  size  of  tubing  available  for  the 
sample  sleeves.  The  inner  diameter  of  the  sample  sleeving  used  had 
to  be  slightly  greater  than  the  driving  face  of  the  transducer.  Once  a 
sleeving  size  was  selected,  the  sample  chamber  size  was  based  upon 
the  outer  diameter  of  the  sleeving.  The  length  of  the  chamber  was 
chosen  to  avoid  finite  wavelength  effects  for  frequencies  as  high  as  one 
kilohertz. 


Figure  17.  Cutaway  View  of  Complete  Coupler. 


The  sample  chamber  and  the  driver  chambers  are  joined 
together  and  kept  aligned  via  alignment  pins.  The  entire  coupler  can 
be  pressurized.  The  pressure  in  the  sample  chamber  is  the  same  as 
the  pressure  in  the  driver  chambers,  minimizing  the  difference  in 
pressure  across  each  driver  diaphragm.  Communication  between  all 
chambers  is  via  a  vent  line.  The  system  is  also  capable  of  being  placed 
under  a  vacuum  to  degas  the  fluid  in  the  sample  and  remove  bubbles 
acquired  while  filling  the  sample  cavity.  O-rings  are  used  to  seal  the 
chamber  against  leakage  of  gas  from  the  driver  side  and  leakage  of 
liquid  from  the  sample  side.  Samples  are  enclosed  in  a  copper  sleeve 
for  easy  installation  and  removal.  The  sleeve  is  three  inches  long  and 
has  an  inner  diameter  of  1.50  inches  and  an  outer  diameter  of  1.66 
inches. 

2.  Fluid  Management  System 

Figure  18  represents  the  fluid  system  used  to  support  coupler 
operation.  A  description  of  fluid  system  operation  is  given  below. 

To  fill  the  chamber,  the  pressure  source  (pressurized  nitrogen 
gas)  is  isolated.  The  chamber  and  the  reservoir  are  then  exposed  to  a 
vacuum.  The  chamber  and  fluid  are  evacuated  to  degas  as  much  as 
possible.  The  fluid  reservoir  is  then  isolated  from  the  vacuum  on  the 
main  vent  side.  The  main  vent  valve  is  then  cracked  open.  The  fluid 
is  driven  by  atmospheric  pressure  into  the  sample  chamber. 

The  vent  line  on  the  chamber  is  constructed  of  high  pressure 
nylon  braided  tubing.  A  fluid  level  is  established  at  the  lower  portion 
of  this  tubing,  just  above  the  connection  to  the  chamber.  The  main 
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vent  valve  is  shut,  and  the  reservoir  is  again  exposed  to  vacuum.  The 
level  in  the  chamber  vent  line  may  lower  out  of  sight,  so  the  fill 
procedure  is  repeated  until  the  level  remains  visible  in  the  tubing. 

Once  a  fluid  level  is  established,  the  chamber  is  isolated  from 
the  vacuum  and  reservoir.  The  valve  to  the  pressure  source  is  opened 
and  pressure  is  slowly  brought  up  to  two  atmospheres.  Pressure  is 
held  there  for  30  minutes.  This  ensures  that  any  remaining  gases  are 
forced  back  into  solution.  The  liquids  to  be  used  in  this  experiment 


Figure  18.  Fluid  System. 


are  Fluorinert  liquids.  These  fluids  are  manufactured  by  the  3M 
Company  [Ref.  16].  This  type  of  fluid  has  been  used  successfully  in 
similar  experiments. 
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IV.  PRELIMINARY  EXPERIMENTS  AND  RESULTS 

A.  INTRODUCTION 

Several  preliminary  experiments  were  preformed  to  validate  the 
concept  that  the  acoustic  load  presented  to  the  face  of  each  driver, 
particularly  a  moving  mass,  can  be  extracted  from  its  input  electrical 
impedance.  Experiments  were  performed  with  a  water-filled  tube 
open  in  air,  with  the  empty  coupler  and  with  the  coupler  filled  with 
water.  One  important  result  of  these  preliminary  experiments  was  the 
need  to  include  a  correction  factor  for  the  difference  in  the  cross- 
sectional  area  of  the  sample  chamber  and  the  driver  diaphragm. 

B.  WATER-FILLED  TUBE 

1.  Description  of  Tube 

Once  we  had  confidence  in  our  ability  to  characterize  the 
driver,  an  experiment  was  conducted  with  a  water-filled  tube  placed 
over  the  driver  to  see  if  the  mass  of  water  added  could  be  extracted 
from  its  input  electrical  impedance.  A  driver  was  fitted  with  a  special 
water-tight  cover  plate  with  tube,  as  shown  in  Figure  19.  The  tube 
and  cover  were  machined  to  the  exact  diameter  of  the  diaphragm, 
including  surround. 

2.  Resonance  Experiment 

A  known  volume  of  water  was  added  to  the  tube,  and  the  shift  in 
resonance  frequency  was  recorded.  Table  4. 1  lists  the  resonance 
frequencies  for  the  volume  of  water  added.  The  lowering  of  the 
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resonance  with  mass  added  is  expected.  These  values  will  be 
quantitatively  analyzed  further  to  determine  if  the  resonance 


Figure  19.  Initial  Test  Apparatus. 

frequency  shift  is  that  expected,  i.e.,  does  the  resonance  frequency 
shift  predict  the  correct  mass  added? 


TABLE  4.1.  RES.  FREQUENCY  FOR  ADDED  WATER. 

Volume  Added  (ml) 

Res.  Frequency  (Hz) 

0.0 

1000 

5.0 

319 

10.0 

226 

15.0 

188 

20.0 

160 

3.  Cross-Sectional  Area  Correction  Factor 

The  resonance  frequency  data  were  reduced  by  comparing  the 
results  of  the  observed  resonance  frequency  with  a  predicted  value.  A 
ratio  of  Equation  3.8  and  3.9  was  established,  and  from  the  known 
stiffness  and  mass  of  the  driver,  the  added  mass  was  determined. 
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When  the  calculated  masses  were  compared  to  the  masses  added,  they 
were  half  of  the  expected  value  for  all  volumes  added. 

The  data  were  analyzed  in  another  manner  by  plotting  the 
period  squared  versus  the  mass  added  and  is  shown  in  Figure  20.  As 
discussed  in  Section  III.B.3.C.,  the  ratio  of  intercept  to  slope  yields 


Correction  Factor  Plot 


Figure  20.  Correction  Factor  Data 


the  mass  of  the  driver.  When  this  was  accomplished,  the  calculated 
moving  mass  was  double  the  known  value. 

This  apparent  discrepancy  in  the  results  can  be  understood  as 
follows.  It  was  assumed  the  driver  acts  as  a  piston  with  uniform 
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displacement  across  the  width  of  the  tube  for  a  specific  current  and 
frequency.  In  reality,  the  edges  of  the  speaker  are  fixed  so  this 
assumption  is  crude  at  best. 

A  better  approximation  of  the  driver  displacement  profile  is 
that  shown  in  Figure  21.  The  displacement  is  assumed  to  increase 


1 
8 
T 

Figure  21.  Assumed  Diaphragm  Displacement. 

linearly  across  the  surround  and  then  stay  constant  across  the  dome. 
Using  this  profile,  the  ratio  of  the  average  velocity  of  the  fluid  above 
the  driver  was  compared  to  the  velocity  of  the  voice  coil,  which  is 
attached  to  the  dome.  The  kinetic  energy  of  the  water  added  was  set 
equal  to  the  kinetic  energy  imparted  by  the  voice  coil.  Referred  to  the 
voice  coil,  the  effective  mass  of  water  added  is  given  by: 


m 


e  ~ 


2  KE 
2  ’ 


(4.2) 


where  KE  is  the  kinetic  energy  of  the  energy  of  the  water,  and  vvc  is 
the  velocity  of  the  voice  coil.  The  ratio  of  me  to  m^ater  depends  only 
on  the  diameters  of  the  tube  and  diaphragm  and  the  diaphragm 
velocity  profile.  We  call  this  ratio  the  "area  correction  factor"  since  it 
is  essentially  the  square  of  the  ratio  of  the  tube  cross-sectional  area  to 
the  cross-sectional  area  of  the  dome.  The  resulting  ratio  for  the 
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velocity  profile  in  Figure  21  is: 

nip 

- —=0.71.  (4.3) 

m  water 

In  other  words,  provided  the  displacement  profile  is  the  one 
assumed,  the  effective  mass  of  the  water  presented  to  the  transducer 
is  only  0.71  times  the  actual  mass  of  water  added.  Since  the  actual 
profile  of  the  transducer  displacement  is  not  really  known,  the 
observed  value  of  0.5  obtained  from  the  resonance  experiments  can  be 
considered  a  "calibration"  value  of  me/m-^ater  for  the  tube.  A  similar 

calibration  will  be  needed  for  the  drivers  in  the  coupler. 

4.  Broad  Band  Experiment 

Once  the  area  correction  factor  for  the  tube  was  obtained, 

experiments  were  conducted  in  an  attempt  to  extract  the  radiation 
impedance  (Zr  )  presented  to  the  diver  over  as  broad  a  frequency 

range  as  possible  (hopefully  100-1000  Hz).  Referring  to  Figure  8,  the 
procedure  is  to  obtain  the  electrical  input  impedance  of  the  driver  and 
the  tube  with  no  added  fluid  (with  single  driver  operation,  the  current 
is  no  longer  halved).  Assuming  negligible  air  loading,  this  can  be 
considered  as  the  series  impedance  of  ZE  +  ZM  .  The  blocked 

electrical  impedance  (ZE)  is  subtracted  from  the  input  electrical 

impedance  and  the  results  are  stored  in  the  HP-4194A.  The  result  is 
the  electrical  equivalent  mechanical  impedance,  ZM  . 

Fluid  is  now  added  to  the  tube  and  the  procedure  repeated. 
The  result  is  the  parallel  combination  of  ZM  and  ZR  .  Since  ZM  is 
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known,  Zr  can  be  solved  for  (Section  II. D.)  using  register 
manipulation  on  the  HP-4194A.  Zr  is  then  obtained  from: 


Z 


r  “ 


(4.4) 


The  real  part  of  Zr  is  the  radiation  resistance,  while  the  imaginary  part 
is  the  radiation  reactance. 

For  the  case  of  fluid  in  the  tube  above  the  voice  coil,  the 
magnitude  of  the  radiation  reactance  is  equal  to  o)Me  where  ©  is  the 
angular  frequency  and  Me  is  the  effective  mass  added.  The  HP-4194A 
has  the  ability  to  perform  the  scalar  multiplication  required  to  extract 
the  added  mass.  This  mass  value  can  be  compared  to  the  mass  added 
to  evaluate  the  effectiveness  of  detecting  the  added  mass  using  the 
acoustically  "soft"  transducers  selected. 

Using  the  area  correction  factor  mentioned  earlier,  the  values  of 
measured  mass  versus  actual  mass  agreed  to  within  five  per  cent  of 
each  other  over  80  per  cent  of  the  frequency  span  between  100  to 
1000  Hz.  Over  the  remaining  20  per  cent  of  the  span  differences 
larger  than  five  per  cent  can  be  attributed  to  the  tube  resonances  that 
occur  in  the  air  in  the  unfilled  portion  of  the  tube.  These  should  be 
eliminated  when  the  chamber  is  tested  full  of  liquid. 
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C.  COUPLER  EXPERIMENTS 

1.  Repeat  of  Water  Pilled  Tube  Experiments  With  Half  of 

Coupler 

Since  the  area  of  the  sample  chamber  is  greater  than  the  area  of 
the  tube  used  earlier,  a  new  area  correction  factor  for  the  drivers 
installed  in  the  coupler  must  be  obtained. 

a.  Resonance  Shift 

The  initial  experiment  conducted  with  the  chamber  was  to 
establish  an  area  correction  factor  for  each  driver,  as  described  in 
Sections  IV.B.2  and  IV. B. 3.  Only  one  driver  is  tested  at  a  time,  with 
the  driver  and  the  chamber  placed  in  the  vertical  position.  Water  is 
added  to  the  open  end  of  the  chamber.  Several  different  volumes  of 
water  are  added  and  the  resonance  frequency  recorded  for  each.  The 
frequencies  and  added  volumes  are  listed  in  Table  4.2.  With  50  ml  of 
water  added,  the  water  level  was  high  enough  to  reach  the  vent  holes 
in  the  sample  sleeving.  Double  this  volume  is  a  good  estimate  of  the 
chamber  volume  with  both  drivers  attached. 


[|  TABLE  4.2.  CHAMBER  CALIBRATION  \\ 

Volume  added  (ml) 

Resonance  Frequency  (Hz) 

0.0 

944.1 

10.0 

266.1 

15.0 

221.3 

20.0 

190.5 

25.0 

171.8 

30.0 

156.7 

35.0 

146.2 

40.0 

138.0 

45.0 

130.3 

50.0 

123.0 

70.0 

104.7 
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b.  Area  Correction  Factor 

The  data  in  Table  4.2  were  reduced  for  the  effective  mass 
added  as  described  in  Section  IV.B.3.  An  area  correction  factor  was 
determined  for  each  driver  based  upon  the  reduced  data.  The  average 
of  the  two  area  correction  factors,  to  be  used  in  further  testing,  is 
0.3108.  The  area  correction  factors  for  the  two  different  drivers 
agreed  to  within  four  per  cent,  reconfirming  the  similarity  of  the  two 
drivers. 

c.  Broad  Band  Experiments 

Once  an  area  correction  factor  was  determined,  broad  band 
measurements  of  input  electrical  impedance  were  made  for  added 
water  volumes  of  10,  20,  and  50  ml.  Applying  the  procedure 
mentioned  in  Section  IV.B.4,  the  data  were  reduced,  resulting  in 
extracted  radiation  mass  as  a  function  of  frequency.  Figure  22  is  an 
example  of  the  results  obtained  for  10  ml  of  added  water.  For  all 
three  added  volumes,  the  extracted  values  of  mass  versus  added 
frequency  were  within  10  per  cent  of  the  actual  value.  Over  most  of 
the  frequency  band,  the  data  was  within  five  per  cent  of  the  correct 
mass. 

2.  Experiments  With  Complete  Coupler 
a.  Air  Filled 

Prior  to  conducting  a  test  of  the  chamber  filled  with  water,  it 
was  first  operated  with  only  air  in  the  chamber.  Both  modes  of 
operation  were  tested:  the  in-phase  "push-push"  mode,  and  the  out- 
of-phase  "push-pull"  mode.  The  resulting  plots  of  impedance 
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magnitude  versus  frequency  are  shown  for  the  two  cases  in  Figures  23 
and  24,  respectively. 


Radiation  Mass  for  10  ml  of  Added  Water 


Figure  22.  Mass  Extracted  from  Impedance  for  10  ml  of  Added  Water. 

An  analysis  of  the  data  reveals  resonances  of  the  air  in  the 
chamber.  The  push-push  data  showed  resonances  at  frequencies  for 
which  the  chamber  length  was  an  integral  multiple  of  the  wavelength. 
The  push-pull  mode  exhibited  resonances  when  the  chamber  length 
was  an  odd  half-integer  multiple  of  the  wavelength.  Also,  these  figures 


demonstrate  the  drivers  do  indeed  appear  electrically  as  one  driver 
with  half  the  impedance,  when  operated  in  parallel.  Since  the  results 


Push-Push  Mode  With  Air  in  Chamber 


Figure  23.  Push-Push  Operation  for  Air  Only. 


were  favorable  for  the  air-only  case,  the  final  experiment  performed 
was  to  fill  the  chamber  with  degassed  water  and  from  the  measured 
input  electrical  impedance  of  the  drivers  wired  in  parallel  out-of¬ 
phase,  extract  the  mass  of  the  water. 
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Push-Pull  Mode  With  Air  in  Chamber 


100  1000  10000 

Frequency  (Hz) 

Figure  24  Push-Pull  Operation  for  Air  Only. 

b.  Water  Filled 

With  the  chamber  filled  with  water,  the  HP-4194A  was  used 
to  measure  the  input  electrical  impedance  (as  a  function  of  frequency) 
of  the  drivers  wired  in  parallel  out-of-phase.  The  acoustic  radiation 
impedance  ZR  was  obtained  as  described  in  Section  IV.B.4.  ZR  was 
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further  reduced  to  obtain  the  radiation  mass  of  the  water  as  described 
in  Section  II. D. 

The  results  are  shown  in  Figure  25.  The  expected  value  of 
radiation  mass  was  50  grams.  The  agreement  is  poor.  At  about  200 
Hz,  there  is  evidence  of  an  unexpected  resonance  in  the  system.  A 
preliminary  analysis  of  the  system  leads  us  to  believe  this  is  a 


"Push-Pull"  Radiation  Mass 


Frequency  (Hz) 


Figure  25.  Radiation  Mass  of  Fluid  Filled  Coupler. 


resonance  of  the  mass  of  fluid  in  the  vent  hole  and  the  spring  of  the 
air  in  the  tubing  above  it.  It  is  expected  that  this  unwanted  resonance 
can  be  eliminated  by  installing  a  valve  in  the  vent  hole  which  can  be 
closed  after  filling  with  fluid. 


56 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  has  been  described  by  which  the  complex  bulk  modulus 
B  and  complex  effective  mass  density  p  for  a  fluid  contained  in  a  rigid 
porous  solid  can  be  extracted  from  measured  input  electrical 
impedance  of  a  pair  of  identical  transducers  wired  in  parallel  in-phase 
and  out-of-phase. 

A  pair  of  moving  coil  transducers  was  selected  and  their 
electromechanical  performance  was  characterized. 

A  cylindrical  brass  coupler  chamber  was  designed  and  built. 

Several  preliminaiy  experiments  were  conducted  to  verify  the 
ability  to  extract  the  density  from  measured  input  electrical 
impedance.  The  results  for  the  most  part  were  accurate  to  within  five 
per  cent  of  expected  values.  The  only  measurements  that  exceeded 
this  limit  were  those  taken  during  preliminary  parallel  driver 
operation  for  the  water  filled  chamber.  Complete  characterization  of 
the  fluid  filled  chamber  still  remains  to  be  done. 

Other  areas  for  future  research  include  adding  a  valve  to  the  vent 
line,  repeating  the  water  filled  experiment,  then  interpret  the  results 
for  both  push-pull  and  push -push  operation.  Once  system  operation 
has  been  verified  with  water,  B  and  p  of  simple  rigid  porous  solids  can 
be  obtained,  providing  additional  validation  of  the  theory  prior  to  the 
testing  of  the  more  complex  porous  solids. 
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APPENDIX 


FLUID  FILLED  POROUS  MEDIA  EFFECTS 


'Hiis  appendix  presents  a  discussion  of  the  effect  of  placing  a  fluid 
in  a  rigid  porous  solid  on  its  density,  bulk  modulus  and  speed  of 
sound.  Derivations  for  a  simple  bulk  fluid  are  followed  by  those  for  a 
fluid-filled  porous  solid.  The  frame  of  the  porous  solid  is  assumed  to 
be  rigid.  The  losses  in  the  fluid  saturated  porous  solid  are  assumed  to 
be  associated  with  the  shear  viscosity  of  the  fluid. 

The  continuity  equation  for  a  bulk  fluid  is  [Ref.  17:pp.  106-108]: 


jjj:  +  v(pu)  =  0. 


(A.  1) 


For  a  fluid  contained  in  a  rigid  porous  solid,  this  formula  becomes 


^  + V«p5>]  =  0, 


(A.  2) 


where. 


(  p)  =  the  spatial  average  of  the  microscopic  fluid  mass 
density. 

(pu)  =  the  spatial  average  of  the  microscopic  fluid 
momentum  density. 

Assuming  a  one  dimensional,  time  harmonic  solution  to  these 
equations,  they  become  in  linearized  form: 


)<0p  +  P0^=°. 


](»<p)+  Po^^  =  °. 


(A.  3) 


(A.  4) 
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Now  consider  Newton's  second  law  for  a  fluid  [Ref.  17:pp.  111-115): 


p^+  (u  *v)  u  =  -Vp  +  r|VxVxu 
at 


(A.  5) 


Assuming  a  one  dimensional,  time  harmonic  solution  and  neglecting 
viscosity  for  the  bulk  fluid,  this  becomes  in  linearized  form, 
)®p°u=-^.  (A.  6) 

The  corresponding  equation  for  a  fluid -saturated  porous  solid, 
including  viscosity  is: 


j®Peff<u>  =  -  Rflow<u>  • 


(A.  7) 


where 

Peflf  =  the  effective  fluid  mass  density, 

<u>  =  the  spatial  average  of  the  microscopic  fluid  velocity, 
<p)  =  the  spatial  average  of  the  microscopic  pressure, 
Rflow  =  the  flow  resistance. 
pe ff  and  Rflow  are  operationally  defined  by: 


^Peff<\^/>  =  kinetic  energy  density  within  the  fluid,  and 

£ 

(u)  Rflow  =  rate  of  decay  of  fluid  kinetic  energy  density 
due  to  shear  viscosity. 

In  general,  peff  £  p.  Also,  for  steady  flow,  Rflow  s  Rpc  and 


Rdc  ^ - — - g’  (A.  7a) 

(pore  size) 

where  r\  is  the  fluid  shear  viscosity. 

It  is  convenient  to  combine  the  velocity  terms  in  Equation  A.7  and 
define  the  complex  fluid  mass  density,  p,  by: 


59 


(A.  8) 


Hence 


P  =  Peff  + 


Rfl 


ow 


(A.  9) 


To  develop  a  wave  equation,  first  the  time  derivative  of  the 
continuity  equation  is  subtracted  from  the  spatial  derivative  of 
Newton's  second  law.  Applying  this  to  the  bulk  fluid  case. 


2  a2p 


to  p  = 


(A.  10) 


ax' 


For  the  porous  solid,  the  corresponding  equation  is: 


2  p  /  \  d  (py 

-co  f-<p>=  — 

Po  ax2 


(A.  11) 


The  next  step  is  to  momentarily  examine  the  bulk  modulus.  The 
definition  of  bulk  modulus  is  [Ref.  18): 

B=Po(||)  (A.  12) 

For  the  fluid  in  the  porous  solid  the  corresponding  definition  is: 

BSp0(*g).  (A.  13) 

Substituting  these  expressions  into  Equations  A.  10  and  A.  11,  the 
following  results: 


2 

-03  p 


_  B  3j> 
‘Po3x2 


for  the  bulk  fluid,  and 


(A.  14) 
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(A.  15) 


-«2<p) 


_ea2» 

p  dx2 


for  the  fluid  saturated  porous  solid. 

Equation  A.  14  is  the  usual  wave  equation  for  the  propagation  of 
sound  in  a  bulk  fluid,  with  a  wave  speed  given  by: 


2  B 

c  =  — . 
P 


(A.  16) 


Equation  A.  15  is  the  wave  equation  for  the  propagation  of  sound  in  a 
fluid-saturated,  rigid  porous  solid,  with  a  wave  speed  given  by: 


(A.  17) 


Hence  the  speed  of  sound  in  a  fluid  saturated,  rigid  porous  solid  is  a 
complex  quantity  which  is  obtained  from  the  complex  bulk  modulus 
and  complex  fluid  mass  densities  in  exactly  the  same  way  the  speed  of 
sound  in  a  bulk  fluid  is  obtained  from  the  bulk  modulus  and  fluid  mass 
density. 
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